Radiation Damping in the Photoionization of Fe^'^"'" 
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A theoretical investigation ol photoabsorption and photoionization of Fe'^*^ extending beyond an 
earlier frame transformation R-matrix implementation is performed using a fully-correlated, Breit- 
Pauli R-matrix formulation including both fine-structure splitting of strongly-bound resonances 
and radiation damping. The radiation damping of 2p — ^ nd resonances gives rise to a resonant 
photoionization cross section that is significantly lower than the total photoabsorption cross section. 
Furthermore, the radiation-damped photoionization cross section is found to be in good agreement 
with recent experimental results once a global shift in energy of « —3.5 eV is applied. These 
findings have important implications. Firstly, the presently available synchrotron experimental data 
are applicable only to photoionization processes and not to photoabsorption; the latter is required in 
opacity calculations. Secondly, our computed cross section, for which the L-shell ionization threshold 
is aligned with the NIST value, shows a series of 2p nd Rydberg resonances that are uniformly 
3-4 eV higher in energy than the corresponding experimental profiles, indicating that the L-shell 
threshold energy values currently recommended by NIST are likely in error. 

PACS numbers: 34.80.Lx,31.15.vj,32.80.Zb 
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Most of what wc know about our universe is obtained 
through spectroscopy. We study either the emission of 
hot plasma sources or the absorption of intervening gas 
between us and a bright object. For example, it is pos- 
sible to probe the muUiphase interstellar medium (ISM) 
through the observation of atomic absorption hnes and 
edges in the spectra of background sources, and the ab- 
sorption properties of the plasma depend on the photoab- 
sorption cross sections of the species present [l| . To inter- 
pret emission spectra, on the other hand, one first needs 
to determine the ionization structure of the gas, which, 
for photoionizcd plasmas, depends on the photoioniza- 
tion cross sections of the chemical elements present. Un- 
hke spectra in the visual region dominated by valence 
electron processes in low ionization species, extreme UV 
and X-ray spectra show inner L- and K-shell processes 
for which photoionization and photoabsorption are signif- 
icantly different. The difference between these two pro- 
cesses arises from the fact that absorption of a photon by 
the atom to an autoionizing state (above the threshold) 
does not necessarily lead to electron emission. 

For the present case of photons incident on Fe^^"*" in the 
vicinity of the L-edge, the difference between photoab- 
sorption and photoionization can be understood by con- 
sidering the prominent 2p — ^ nd absorption resonances. 
Following absorption of a photon. 
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an intermediate resonance state can either decay via au- 
toionization, or Auger decay, leading to ionization of 
Fe"'^^+ into the Fe^^"*" ion plus a free electron, or it can 
radiatively stabilize, leading to no charge change of the 



2p^3s^nd 



2/3s + e- 
2p^nd + e: 
2/3s2 + hv' (AP° 



(A^P'^"* ind. of n) 
-n-3) 



(2) 
(3) 
(4) 



— > 2p^'isnd+hv' (Af'^^* ind. of n) (5) 

Thus, the alternative decay pathways of radiative stabi- 
lization redirect some of the initial photoabsorption am- 
plitude, giving a reduced, or damped photoionization 
cross section relative to the photoabsorption cross sec- 
tion. 

The extent of this so-called radiation damping effect, 
for an isolated resonance, is essentially given by the au- 
toionization branching ratio 



0"PJ 
CPA 



(6) 



where upi and a pa ai'c the photoionization and pho- 
toabsorption cross sections, respectively. The total au- 
toionization and radiative rates, A*°* and v4*°*, are each 
the sum of contributions for which the nd Rydberg elec- 
tron participates, with partial rates scaling as n~^ , and 
for which the Rydberg electron is a spectator, with par- 
tial rates that arc independent of n. Thus, the spectator 
Auger and radiative decay branches (Eqs. [3] and [5]) dom- 
inate as n — > oo. On the other hand, due to angular 
momentum, or geometric, and radial, or dynamic, con- 
siderations of the 3d orbital compared to the 3s orbital, 
the 2p^?)S^2>d resonance decays predominantly via partic- 
ipator pathways, as given in Eqs. [2] and SI 

To date, essentially all of the calculated photon-atom 
data for multiply-charged ions has been produced by the- 
oretical photoionization calculations without considering 
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the alternate radiative decay pathway in Eqs. |4] and [5] 
(see, for example, Refs. [1, Q)- These undamped cross 
sections are equivalent to the photoabsorption cross sec- 
tions regarding absorption strength (although the reso- 
nance width is underestimated), and thus it had always 
been assumed that api = crpA- For Fe^^^, perhaps the 
most definitive calculations to date are the R-matrix pho- 
toionization calculations in both LS coupling [5| and us- 
ing an LS-coupled to a JK-coupled frame transforma- 
tion (LSJKFT) method, along with resonance broad- 
ening in order to correctly characterize the resonances 
0. In that latter study (hereafter referred to as Pa- 
per I), it was demonstrated that the LSJKFT method 
was capable of reproducing all fine-structure splitting 
effects of a full Breit-Pauli R-matrix (BPRM) calcula- 
tion except for the lowest 2s^2p^3s^3(i and 2s2p^3.s^3p 
resonances. This is because in the LSJKFT method, 
an extremely efficient computational approach compared 
to the BPRM method, the multi-channel quantum de- 
fect (MQDT) equations for the outer-region solutions 
are modified, thereby incorporating fine-structure split- 
ting only for those resonances that reside outside the 
R-matrix " box" ; those that reside within the R-matrix 
region (i.e., all states described by orbitals only up to 
n = 3, such as the 2s^2p^3s^3d resonances) did not in- 
clude fine-structure splitting effects. 

The agreement between LSJKFT and BPRM meth- 
ods (except for the strongly-bound resonances) was 
first demonstrated only for a minimal configuration- 
interaction (CI) case, since the large-scale CI calcula- 
tion which would have been required for a more con- 
verged calculation and which was performed only within 
an LSJKFT approach, would have been prohibitively 
large to undertake within a full BPRM approach. Inter- 
estingly, it was found that, in terms of general qualitative 
features, the large-CI LSJKFT cross section, convoluted 
with a sufficiently broad energy distribution, was essen- 
tially the same as that from a minimal-CI LS calculation. 
Thus, it is necessary to include higher-order CI and fine- 
structure effects only if a detailed resonance description is 
desired. It should be noted that all of those earlier calcu- 
lations were for the total photoabsorption cross section, 
and did not consider radiation damping effects. Further- 
more, those earlier cross sections were all preconvolutcd 
with a constant Lorcntzian width of 0.1 Ryd Q, so 
a detailed investigation of the resonance profiles was not 
performed. 

More recently, an electron-beam ion trap (EBIT) was 
used at a synchrotron facility to detect the final pro- 
duction of Fc^'^"*" ions thereby measuring the pho- 
toionization cross section at an energy resolution as low 
as 150 meV (see also Ref. [l0|)- Thus, it is now desir- 
able to compute radiation- damped photoionization cross 
sections, since only total photoabsorption cross sections 
were reported in Paper I. Furthermore, to compare di- 
rectly with the finer resolution of the experimental data, 
it is necessary to consider higher-order CI effects and 
fine-structure splitting for all resonances, including the 



lowest 2s^2p^is'^?,d resonances for which a frame trans- 
formation approach is inapplicable. The present study 
improves upon the theoretical work of Paper I and other 
earlier calculations [ll| by using the same large-scale 
CI description of Paper I within a full BPRM calcula- 
tion, which treats fine-structure splitting correctly for 
even the lowest resonances, and by including radiation 
damping effects [12|. The BPRM calculations arc per- 
formed using the standard UCL/Belfast R-matrix codes 
with additional modifications for incorporating radia- 
tion damping (l^. as described in earlier electron- impact 
excitation [isj and dielectronic recombination (l3 | stud- 
ies. A brief description of the present methodology is 
given below, followed by a detailed comparison between 
these new results and the recent experimental results. 

The atomic structure used is essentially the same as 
that of Paper I. The Fe^^"*" target states are described by 
an orbital basis set that is generated from a Hartrec-Fock 
[isj calculation on the ls^2s^2p^3s ground state configu- 
ration, with additional 3p and Zd orbitals generated from 
frozen-core Hartree-Fock calculations on the ls^2s^2p^3p 
and ls^2s^2p^id excited configurations, respectively. All 
target states of the form 2s^2p^3Z, 2s^2p^3s3/, and 
2s2p^3s3Z are used in the close-coupling expansion and 
are described by a larger CI basis consisting of all con- 
figurations consistent with single and double promotions 
out of those three base configurations. 

The corresponding (A'^-l-l)-electron bound, continuum, 
and resonance states of Fe^'*+ are described by a basis 
consisting of single, double, and triple promotions out 
of the 2s^2p^3s^ base configuration. Lastly, unlike in 
Paper I, corrections to the A^-electron Hamiltonian are 
applied in order to realign the target energies with the 
recommended NIST [l^ values (see Table I), and the 
binding energy of the Fe^^+ ground state, relative to the 
Fe^^"*" state, is also aligned with the NIST value. Thus, 
of particular importance, the present theoretical photon 
energy thresholds which we report are aligned with the 
corresponding NIST values. 

The important details regarding radiation damping are 
that ah 2s22/3s2 2s'^2p^2,sM (J = 0,1,2) bound 
states are included as final type II (inner-region) radia- 
tive decay states, resulting in a complex inner- region R- 
matrix jl3 - [l4l | . This accounts for all the participator ra- 
diative decay pathways in Eq. 2] and spectator radiative 
decay of the 2^^3523^ resonance. The spectator radiative 
decay of all higher nd resonances in Eq. [SJ on the other 
hand, is accounted for by a modification to the MQDT 
equations in the outer region [l^ - [l^ . where the type I 
(3s — ;> 2p core) decay width used in the E ^ E — iTrl2 
modification is = 0.45 meV. 

The target states included explicitly in our R-matrix 
calculation only account for the participator Auger decay 
of Eq. [5] and the n = 3 spectator Auger decay of Eq. [3] 
to the so-called main line photoionization continua. The 
satellite photoionization continua of Eq. [3] for n > 3 are 
instead included in our formulation via an optical poten- 
tial MQDT modification procedure [l^, similar to that 
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done for the spectator radiative decay. This presents a 
difficulty in extracting a total theoretical photoionization 
cross section, since the above described methodology for 
including radiation damping and n > 3 spectator Auger 
decay yields only total photoabsorption and main-line 
photoionization cross sections The difference be- 

tween the two includes both satellite photoionization and 
radiative damping amplitudes. As ri — ?> oo, the branch- 
ing of the two can be extracted by using the spectator 
rates in Eq. [HI whereas for low-n, the difference is due 
purely to radiation damping. For intermediate n, how- 
ever, the branching requires a deeper investigation be- 
yond the scope of the present work. 

In Fig. 1, the present BPRM results are compared to 
the earlier results of Paper I and to recent experiment Q . 
In Fig. la, our new BPRM photoabsorption results, con- 
voluted with a Lorentzian width of 0.1 Ryd, are found 
to be qualitatively similar to the earlier LSJKFT re- 
sults, but differ quantitatively in three minor respects. 
First, there is a uniform energy difference between the 
two since the earlier results were computed using the- 
oretical thresholds (NIST values were not available for 
the 2p^3s^ autoionizing states at that time) whereas our 
present results have been aligned to the NIST thresh- 
olds. Second, the LSJKFT results were preconvolutcd 
0] with an 0.1 Ryd Lorentzian Q only above 931.55 
eV, whereas the present BPRM results have preconvo- 
lutcd even the lowest resonances. Third, and most im- 
portantly, the LSJKFT results do not account for fine- 
structure splitting of the lowest 3d resonances, as noted 
earlier. The BPRM and LSJKFT cross sections are also 
compared to each other in Fig. lb on a linear scale, con- 
voluted with a Voigt profile (Lorentzian and Gaussian 
widths Tl = 2.11 cV and Tq = 4.0 eV, respectively). 
It is seen that, except for the slight global energy shift 
and the BPRM fine-structure splitting of the 25^2^53523^ 
resonances, the two cross sections are nearly identical. 
Thus, the only improvement over the LSJKFT method 
by using instead a BPRM method is the inclusion of fine- 
structure splitting of the lowest resonances, as had been 
noted in Paper I. 

Using the same resolution, we compare our total pho- 
toabsorption cross section to the computed main-line, 
damped photoionization cross section in Fig. Ic. As 
noted above, the main-line photoionization cross section 
for the n = 3 resonances below 950 eV photon energies, 
and, for the most part, the n = 4 resonances, accounts 
for the total photoionization cross section, and thus there 
is appreciable 3d — ^ 2p radiation damping of the lowest 
2p^3s^nd resonances; the 2s2p^3s^3p resonance at 880 
eV, on the other hand, is not damped significantly due 
to the much smaller 3p — t- 2s radiative rate. 

We compare the main-line photoionization cross sec- 
tions to the experimental results (arbitrarily scaled) in 
Fig. Id. Except for an approximately constant energy dif- 
ference, good agreement between the two results is found 
in resonance profiles up to about 1050 cV. It is interesting 
to note that the relative heights of the 2p^3s^ P::^/2)3d 
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FIG. 1: (Color online) Fe^^+ photoabsorption and main-line 
photoionization cross sections in the vicinity of the L-edge. a) 
Earlier LSJKFT photoabsorption results [y] (blue line) and 
present BPRM photoabsorption results (green line), using a 
Lorentzian width of 0.1 Ryd. b) Same as a) using a Voigt 
profile (Fl = 2.11 eV, To = 4.0 eV). c) Present BPRM pho- 
toabsorption results (green line) and main-line photoioniza- 
tion results (red line) using the same Voigt profile, b) Present 
main-line photoionization results (red line) and experimental 
results (black data points). 

and 2p^3s^{^Pif2)'^d resonances are the same in both the 
BPRM photoionization and experimental cross sections, 
whereas the BPRM photoabsorption cross section indi- 
cates a much different height ratio. 

The comparison of relative 2p^3s^3d resonance heights 
would seem to validate experimentally the extent of 
radiation damping. Nevertheless, we can get a more 
quantitative assessment by examining absolute reso- 
nance strengths. By fitting our photoionization and 
photoabsorption cross sections, in the vicinity of the 
2p53s^(^Pi/2)3rf resonance, to an energy-normalized 
Lorentzian multiplied by a strength factor, we obtained 
integrated resonance strengths of 130 Mb-eV and 200 
Mb-eV, respectively. Table I indicates that our pho- 
toionization resonance strength is indeed in agreement 
with the experimental value of 110±60 Mb-eV whereas 
the photoabsorption resonance strength is outside of the 
experimental uncertainty, providing further validation of 
the extent of radiation damping. 
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TABLE I: 2p^3s^(^Pi/2)3d Resonance Strengths (Mb-eV). 



Experiment ^ 


110 i 60 


R-matrix Photoionization 


130 


R-matrix Photoabsorption 


200 


TABLE IL Fe"+ Resonance Enerj 


;ies and Fe^^+ Thresholds 


(eV). 





Level 


Exp. [9] 


Present 


NIST fl6l 


Exp. [18] 


Fe^*+(2p='3s^(^P3/2)3d) 
Fe"+(2p=3s2(^P3/2)4d) 
Fe"+(2p=3s2(^P3/2)5d) 


794.7 
964.3 
1040.9 


798.77 
967.81 
1044.22 






Fei=+(2p=3^^(^P3/2)) 




1172.76° 


1172.76 


1170.88 


Fe^''+(2p''3s^(^Pi/2)3d) 
Fe"+(2p=3s^(^Pi/2)4d) 
Fe"+(2p=3s"( = Pi/2)5d) 


807.1 
976.0 
1053.5 


811.13 
979.25 
1055.9 






Fei=+(2p=3s=^(^Pi/2)) 




1185.16° 


1185.16 





"Theoretical thresholds aligned to NIST values [l^ . 



The comparison m Fig. Id also indicates that the 
BPRM resonance energies are uniformly higher than 
the corresponding experimental values, even though the 
BPRM cross sections have been aligned to the n ^ oo 
thresholds of the 2p^3s'^(^P3/2)nd and 2p^3s'^(^Pi/2)nd 
Rydberg series given by NIST as 1172.77 and 1185.17 
eV, respectively. Table II lists the resonance energies, 
indicating that the resonance energy difference does not 



approach zero with the expected behavior, assum- 
ing the same n — >■ cx) Rydberg limits. This seems to 
suggest that our threshold values used - those recom- 
mended by NIST [igl. [20j - are too high. Indeed, a more 
recent experiment [l8| and subsequent theoretical analy- 
sis [U reports a lower 2p^3s'^{'^P3/2) threshold of 1170.88 
eV, which would account for about 2 eV of the differ- 
ence in resonance energy positions. Although newer mea- 
surements for the 2p^3s^(^Pi/2) threshold energy are not 
available, it is likely that this NIST value is also in error. 

In conclusion, we have reported new R-matrix calcula- 
tions that improved on our earlier work of Paper I by in- 
cluding fine-structure effects for the lowest resonance and 
radiation damping effects for all resonances. We find that 
radiation damping is significant for the 2p^3s^(^i-i/2)3(i 
resonance, in particular. Our radiation-damped pho- 
toionization resonance strengths show good agreement 
with the experimental results, but there is a nearly uni- 
form difference in energy positions between the present 
theoretical and the experimental Q resonance positions. 
Since the present calculations have been aligned to the 
NIST experimental threshold value [lO], there is an ap- 
parent inconsistency between theoretical and experimen- 
tal resonance positions which is not fully resolved from a 
newer threshold measurement [l^. 
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